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The sequence of ground states for SrTiOa film subjected to epitaxial strain as well as to me¬ 
chanical stress along the [001] and [110] axes is calculated from first principles within the density 
functional theory. Under the fixed-strain boundary conditions, an increase in the lattice parameter 
of a substrate results in the 14cm —>■ 74/mcm —>■ Ima2 —> Cm —>■ Fmm2 —>■ /ma2(II) sequence 

of ground states. Under the fixed-stress boundary conditions, the phase sequence is different and 
depends on how the stress is applied. It is revealed that the simultaneous presence of competing 
ferroelectric and antiferrodistortive instabilities in SrTiOa gives rise to the appearance of metastable 
phases, whose number increases dramatically under the fixed-stress conditions. In the metastable 
phases, the octahedral rotation patterns are shown to differ substantially from those in the ground 
state. It is suggested that in systems with competing instabilities, each polar phase has its optimal 
octahedral rotation pattern which stabilizes this phase and creates a potential barrier preventing 
this phase to be transformed into other structures. 

PACS numbers: 61.50.Ah, 61.50.Ks, 64.60.My, 77.80.bn, 77.84.Cg 


I. INTRODUCTION 

Elastic strain is widely used today to improve the prop¬ 
erties of electronic materials. For example, the effect of 
the strain-induced decrease of the acceptor binding en¬ 
ergy in germanium enabled to create unique photode¬ 
tectors working in the far infrared.^ An increase of the 
electron mobility in strained silicon enabled to consider¬ 
ably improve the performance of silicon field-effect tran¬ 
sistors,^ and the use of highly strained layers in pseudo- 
morphic high electron mobility transistors made it pos¬ 
sible to deliberately tune the energy diagrams of these 
heterostructures in order to significantly improve their 
char act erist ics 

The appearance of a strain-induced ferroelectricity in 
different dielectrics significantly extends the functional¬ 
ity of these materials and enables to offer new, previ¬ 
ously known designs of electronic devices. The strain 
of incipient ferroelectrics leads to particularly impressive 
results.®^ For example, the stretching of thin films of 
strontium titanate SrTiOa grown on DyScOa substrates 
increases the Curie temperature in this material, which 
is nonpolar in the absence of strain, to ~300 K,IS1 and the 
compression of these films grown on silicon substrates 
increases this temperature up to ~410 K (Ref. [S|. 

Strontium titanate is an incipient ferroelectric exhibit¬ 
ing a competition between the ferroelectric and antifer¬ 
rodistortive (octahedral rotational) instabilities.^ The 
latter of them is the cause of the phase transition ob¬ 
served in SrTiOa at about 105 K. The first work in 
which the influence of strain on the ferroelectric prop¬ 
erties of strontium titanate was considered within the 
phenomenological approach was the work of Uwe and 
Sakudo.E^Pertsev et aZ.l^have developed this approach to 
describe the effect of the epitaxial strain on SrTiOa thin 
films. By expanding the thermodynamic potential in a 
power series of two order parameters (polarization and 
octahedral rotation) to the fourth order, the authors ob¬ 


tained a rich pressure-temperature phase diagram with a 
large number of different phases. For the set of material 
constants of SrTiOa used by the authors, all obtained so¬ 
lutions had the order parameters directed along the axes 
of the cubic structure. Subsequent first-principles cal¬ 
culations in which the antiferrodistortive instability was 
neglect ed,li2HII] however, have shown that in stretched 
SrTiOs films the polarization should be directed along 
the [1101 axis, in agreement with the experimental 
data.™ A more thorough first-principles study of the 
low-temperature phases of SrTiOa, in which both fer¬ 
roelectric and antiferrodistortive instabilities were taken 
into account,^ showed that in highly stretched SrTiOa 
films the polarization is indeed directed along the [110] 
axis, but at low strain a phase with the [100] polar¬ 
ization appears. This orientation of the polarization 
was observed later in anisotropically straine d SrT iO.-j 
films grown on GdScOa and DyScOa substrates. O^ I ^^ I The 
phase diagrams of strained SrTiOa for other sets of mate¬ 
rial constants were studied within the phenomenological 
and the phase-field approaches in Refs. [niETMl 

When studying the strain effects on the properties of 
different materials, the fixed-strain boundary conditions 
are usually used. However, the study of the fixed-stress 
boundary conditions is also important for real systems. 
To explain this, let’s consider a thin epitaxial film which 
has domain structure and is fixed on a substrate. If the 
substrate is incompressible, then the in-plane lattice pa¬ 
rameters in these domains would be equal to the lattice 
parameter of the substrate. However, in real systems 
the lattice parameters in the domains will be different 
because the substrate is compressible, and the difference 
between them will increase with increasing film thickness. 
In the limit of highly compressible (flexible) substrate, 
if one neglects mechanical and electrical boundary con¬ 
ditions at the domain walls, we come to the fixed-stress 
boundary conditions. In practice, the fixed-strain bound¬ 
ary conditions can be realized in ultrathin epitaxial films, 
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whereas the fixed biaxial stress boundary conditions can 
be realized in films grown on flexible substrates or thin 
plates which are fixed at their edges on a set of piezoelec¬ 
tric actuators. The actual boundary conditions in real 
systems are intermediate between the two limiting cases, 
and this is why the analysis of both cases is necessary 
when studying the strain effects. 

In this work, we consider the ferroelectric and antifer- 
rodistortive instabilities in SrTiOa under the fixed-stress 
boundary conditions and compare them with the results 
obtained for the fixed-strain boundary conditions. Two 
different ways of applying stress to the film are exam¬ 
ined: a uniaxial stress normal to the film plane and a 
biaxial stress in the film plane. We show that these two 
ways of applying stress result in different phase diagrams. 
We demonstrate that the simultaneous presence of com¬ 
peting ferroelectric and antiferrodistortive instabilities in 
SrTiOa gives rise to the appearance of previously un¬ 
known metastable phases, whose number increases dra¬ 
matically under the fixed-stress conditions, and explain 
the origin of the metastability. 


II. CALCULATION TECHNIQUE 

Elastic deformation of SrTiOa films was realized in two 
ways. To create the fixed-strain boundary conditions, 
the ( 001 )-oriented film was grown on a cubic substrate 
with the lattice parameter oq, which was varied within 
± 2 % of the lattice parameter of cubic strontium titanate. 
To create the fixed-stress boundary conditions, the stress 
was applied along the 2 ; axis normal to the film plane 
(the pressure is P[ooi]) or in the xy plane (the pressure is 
P[iio]). The film was considered free to relax in directions 
normal to the applied stress or strain. 

Calculations of the equilibrium lattice parameters and 
atomic positions in films were performed within the first- 
principles density functional theory using the ABINIT 
software. The exchange-correlation interaction was de¬ 
scribed in the local density approximation (LDA). The 
pseudopotentials of atoms constructed using the RKKJ 
schem^^ were borrowed from Ref. [5^ The plane-wave 
cutoff energy was 30 Ha (816 eV), the accuracy of self- 
consistent energy calculations was better than 10“^° Ha. 
For the integration over the Brillouin zone, a 8 x 8 x 8 
Monkhorst-Pack mesh for the cubic cell or meshes with 
equivalent density of fc-points for low-symmetry phases 
were used. The relaxation of atomic positions and lattice 
parameters was performed until the forces acting on the 
atoms become less than 2 • 10“® Ha/Bohr (0.1 meV/A). 
Near the boundaries between the phases, the accuracy of 
the relaxation was increased to 2 • 10“^ Ha/Bohr. 

As the energy difference between some phases in this 
work can be as small as 0.1 meV, the accuracy and con¬ 
vergence of calculations are very important issues. Insuf¬ 
ficient /c-point density, parallel shift of all atoms in the 
unit cell, and different sets of A:-points for different struc¬ 
tures are the main sources of errors. The convergence 


studies have shown that an increase of the fc-point den¬ 
sity by a factor of three changes the energies of phases by 
no more than 0.02 meV. More important errors may re¬ 
sult from a parallel shift of all atoms in the unit cell. The 
standard deviation of the energy of phases calculated for 
20 different random parallel shifts was equal to 0.028 meV 
and was independent of the fc-point density. As to the 
third source of errors, in most calculations we used the 
same sets of the fc-points when performing the integra¬ 
tion over the Brillouin zone, but the tests have shown 
that the changes in the energy of phases calculated for 
different sets of fc-points do not exceed 0.03 meV. The er¬ 
rors in the enthalpy calculations resulting from variations 
of the unit cell volume and local stresses were less than 
0.003 meV. So, the accuracy of our calculations enables 
to reliably distinguish the energy difference of 0.1 meV. 

In searching for the ground state, after the relaxation 
of a structure of each relevant phase for each value of ap¬ 
plied stress or strain, the phonon spectrum and the elas¬ 
tic tensor were calculated and it was checked whether all 
optical phonon frequencies at all high-symmetry points 
of the Brillouin zone are positive and whether the de¬ 
terminant and all leading principal minors constructed 
from components of the elastic tensor are positive (the 
stability criterion). If the criterion is not satisfied, small 
distortions corresponding to the least stable phonon were 
added to the structure, and the search for the ground 
state was continued. The technique of phonon spectra 
and elastic tensor calculations was similar to that de¬ 
scribed in Ref. |25l We considered only monodomain 
states because domain walls usually have a positive en¬ 
ergy and so the multidomain solutions have a higher en¬ 
ergy. 

III. RESULTS 

The phonon spectrum calculations for the high- 
temperature Pm3m phase of SrTiOa confirmed a well- 
known result that there are three types of instabilities 
in the phonon spectrum: the ferroelectric one associated 
with the Tis mode at the center of the Brillouin zone 
and two antiferrodistortive instabilities with respect to 
rotations of the oxygen octahedra described by R 25 and 
M 3 modes at the boundary of the Brillouin zone.^^IllSI 
Since the M 3 instability is weak and disappears when 
the mode at the R point is condensed, we will neglect 
it when searching for the ground state. The calculations 
showed that in the absence of strain, the ground-state 
structure of strontium titanate is the /4/mcm phase.^^ 

A. Fixed-strain boundary conditions 

To better understand the influence of the fixed-stress 
boundary conditions on the phase diagrams, we first con¬ 
sider the phase diagrams for SrTiOs films grown on a 
cubic substrate (the fixed-strain boundary conditions). 
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FIG. 1. (Color online) Energies of different phases for SrTiOa 
film grown on a cubic substrate with the lattice parameter 
ao- The energy of high-symmetry PAjirimm phase is taken 
as the energy reference. The vertical line indicates the lattice 
parameter of cubic SrTiOa. Short arrows show the lattice pa¬ 
rameters at which the frequencies of soft optical modes vanish. 


Our calculations (Fig. show that the sequence of the 
ground states [/4cm —>■ lA/mcm —>■ Ima2 Cm ^ 
Fmm2{\l) —>■ /ma2(II)] in this case is unexpectedly com¬ 
plex. The transition between the lAcm and lA/mcm 
phases occurs at ag ^ 7.279 Bohr and is not shown in 
the figure (see Table A1 in the Appendix). At oq « 
7.316 Bohr, the /d/mcm phase becomes unstable and 
transforms to the Ima2 phase, in which the octahedra 
are rotated around the [001] pseudocubic axis, whereas 
the polarization is along the [110] one (Table [I]). For all 
strains, this phase has a lower energy as compared to the 
Fmm2 phase in which the rotations are around the same 
axis, but the polarization is along the [100] axis. The 
transition between the Ima2 and Cm phases occurs at 
ao ~ 7.345 Bohr and is a first-order transition (the molar 
volumes of the two phases differ by 0.75% at the tran¬ 
sition point). Near the phase transition point, the Cm 
phase can be regarded as a slightly distorted Jma2(II) 
phase, ^ but the polarization vector quickly rotates in 
the plane from the [110] to [100] direction with increasing 
ao, and already at ag = 7.38 Bohr the Fmm2(II) phase 
becomes the ground state. At ag = 7.39-7.40 Bohr, the 
energy of the Fmm2(II) phase turns out to be slightly 
lower than that of the /ma2(II) phase (by approximately 
0.1 meV, see Table A1 in the Appendix for the ener¬ 
gies of the phases). It is interesting that both phases 
with the lowest energy [Fmm2(II) and /ma2(II)] in this 
region satisfy the stability criterion (i.e., the former of 
them is stable and the latter is metastable), so that in 
this region the phases with polarizations along the [110] 
and [100] pseudocubic axes can coexist. This, in partic¬ 
ular, enables to explain why two different orientations of 
polar izatio n have been observed in two different experi- 
mentJi^lSI on SrTiOa films grown on DyScOa substrates 


TABLE I. Space groups and order parameters for phases con¬ 
sidered in this work. The first three symbols in the order pa¬ 
rameter denote Cartesian coordinates of the polarization vec¬ 
tor and the last three symbols denote the rotations around 
the pseudocubic axes. The lattice parameters and atomic 
positions for typical representatives of the structures can be 
found in Tables A3 to A12 in the Appendix. 


Space group 

Order parameter 

74/mcm 

(OOOOOd-) 

/4cm 

(OOpOOfli) 

Fmm2 

(pOOOO.))), (OpOOOc)-) 

Ima2 

(ppO OO^i) 

Fmmm 

(000 000), (000 000) 

Fmm2(\\) 

(pOOO0O), (OpO0OO) 

Ima2{ll) 

(ppO 000) 

Cm 

(P 1 P 2 O 01020) 

Fmm2{lll) 

(OOpO0O), (OOp0OO) 

Pc 

(PlPlP3010103) 


and to understand the sensitivity of polarization to the 
anisotropic strain of the substrate. At ag > 7.40 Bohr, 
the /ma2(II) phase becomes the ground state. 


B. Fixed-stress boundary conditions 

Under the fixed-stress boundary conditions, the se¬ 
quence of the ground states resulting from compression 
or stretching of the film is different for two ways of 
applying stress. This follows from the difference be¬ 
tween the equations used to calculate the enthalpy. At 
p II [001], the enthalpy is calculated using the formula 
H = Etot +P[ooi]^e 3 , where /?tot is the total energy, 
P[ooi] is the pressure, V is the unit cell volume, and 
£3 is the strain tensor component (in Voigt notation) 
normal to the film. At p || [110], the enthalpy is 
H = Etot +P[iig]V{ex + € 2 ), where ei and £2 are the 
strain tensor components in the film plane. 

When the stress is applied, the symmetry of the 
high-temperature Pm3m phase is lowered to tetragonal 
PA/mmm phase. As for the low-temperature lA/mcm 
phase, the formation of two structures differing by the 
direction of octahedral rotation axis is possible under 
stress. Because of a noticeable spontaneous strain ac¬ 
companying the octahedral rotations in the unit cell, the 
resulting lA/mcm phase has a lower enthalpy for biaxi- 
ally compressed films and the Fmmm phase has a lower 
enthalpy for biaxially stretched films (Fig. [^. In this 
work, a positive sign of the applied stress means the com¬ 
pression and a negative one means the tension. Since 
the compression of the film along the z axis results in 
its expansion in the xy plane, to facilitate the compar¬ 
ison of results for two ways of applying stress, we used 
two opposite directions of the horizontal axis in Fig. [^ 
The in-plane strain of the film at a maximum pressure 
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FIG. 2. (Color online) Enthalpies of different phases for 
SrTiOa film subjected to (a) uniaxial compression along the 
2 axis and (b) biaxial stretching in the xy plane. The en¬ 
thalpy of high-symmetry P^jmmm phase is taken as the en¬ 
ergy reference. Short arrows indicate the stresses at which 
the frequencies of soft phonons vanish. Solid lines connect¬ 
ing the points refer to stable and metastable phases, whereas 
dashed lines indicate the nnstable phases. Long vertical lines 
with arrows show the directions of structnral relaxations at 
pressnres where the metastable phases become unstable. 


(IpI = 44.1 kbar) is 2.13% for p || [001] and 1.85% for 

pII [110]. 

The existence region of the lA/mcm phase on the phase 
diagram depicted in Fig. [^is limited by the softening of 
the ferroelectric A 2 u mode (under the in-plane compres¬ 
sion) and of the ferroelectric mode (under the in-plane 
stretching). When crossing the boundaries of the exis¬ 
tence region (shown by short vertical arrows in Fig. |^, 
the frequencies of corresponding modes become imagi¬ 
nary, and the symmetry of the unit cell is lowered, re¬ 
spectively, to /4cm and /ma2.1^The same applies to the 
Fmmm phase, in which the softening of the Biu mode 
upon the in-plane compression results in the appearance 
of the Fmm2(III) phase polarized along the 2 axis, and 
the softening of the B^u mode upon the in-plane stretch¬ 


ing results in the appearance of the Fmm2{ll) phase po¬ 
larized along the x axis (we assume that in both cases 
the octahedral rotations are around the y axis). We note 
that for two ways of applying stress we are discussing, the 
values at which the phase transitions occur are very dif¬ 
ferent. For example, for the /4/mcm-/4cm phase tran¬ 
sition the transition pressures are P[ooi] = —6.5 kbar and 
P[iw] = +13.1 kbar. 

While there is only one li/mcm-IAcm phase transi¬ 
tion in the in-plane compressed SrTiOa films, the se¬ 
quence of the ground states in stretched films is much 
more complex. The calculations of the phonon spectra 
of the Fmmm phase show that with increasing in-plane 
stretching not only the polar B^u mode is softened, but 
also the polar B 2 u mode is softened. At a pressure of 
Pliio] = —4.9 kbar or p[ooi] = +18.2 kbar, the former 
mode results in the transition to the Fmm2{lV) phase 
polarized along the x axis. With further increase of pres¬ 
sure (at —12.7 kbar or -1-30.3 kbar, respectively), the 
B 2 u mode induces the transition to the Cm phase, in 
which the polarization is rotated in the xy plane. For 
the in-plane stretching, the Cm phase transforms to the 
/ma2(II) phase at ^30 kbar [see Figj^b) and enthalpies 
of different phases in Table A2 in the Appendix]. In the 
case of the out-of-plane compression, the Cm phase re¬ 
mains the ground state up to at least 45 kbar [Fig. [^a)]. 

The pressure dependence of two components of spon¬ 
taneous polarization for both ways of applying stress 
is shown in Fig. The changes in polarization upon 
the Fmmm —>■ F'mm2(II) and F'mm2(II) —> Cm phase 
transitions and the fact that the polarization abruptly 
changes from the Px Py state in the Cm phase to 
the Px = Py state characteristic of the /ma2(II) phase 
are clearly seen. From the comparison of the polariza¬ 
tions for two ways of applying stress it follows that the 
polarization is about 2.5 times higher for biaxial stretch¬ 
ing than for uniaxial compression, whereas the phase se¬ 
quence is the same. The reason for this is simple: the 
/'[iio]^(£i + £ 2 ) term in the enthalpy makes larger in¬ 
plane strains more favorable under p[iio] < 0. As the in¬ 
plane polarization strongly increases with increasing ten¬ 
sile strain (£ 1 + 62 ), the polarization in biaxially stretched 
films is considerably higher than that in uniaxially com¬ 
pressed hlms. Under the p II [001] compression, the polar¬ 
ization in Fmm2(ll) and Cm phases is much lower than 
that in the metastable /mo2(II) phase because the rota¬ 
tion around the [ 100 ] axis strongly reduces the (ei + € 2 ) 
values. 


C. Hysteresis and metastability effects 

An analysis of the stress dependence of the soft-mode 
frequencies in the ground states of SrTiOa for the fixed- 
stress boundary conditions (Fig. |4) reveals a number of 
previously known features. First, the pressures, at which 
the extrapolated squares of the soft-mode frequencies 
vanish on both sides of the lA/mcm-IAcm phase tran- 
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FIG. 3. (Color online) Components of polarization as a func¬ 
tion of pressure for two ways of applying stress. Data for ten¬ 
sile p II [110] stress are shown by filled symbols and data for 
compressive p || [001] stress are shown by open symbols. The 
dotted and dashed lines show the polarization in the 7ma2(II) 
phase obtained for p || [110] and p || ]001], respectively. 



FIG. 4. (Color online) Frequencies of low-energy phonons in 
the ground-state structures as a function of the in-plane stress 
{p II [110]). Dotted lines show the behavior of frequencies in 
the metastable /ma2(II) phase that coexists with the ground- 
state Cm structure. 


sition, differ by 0.76 kbar for p || [110] (hysteresis loops 
are not shown in the figure). For p \\ [001], the difference 
between these pressures is more significant, 1.88 kbar. 
This means that in both cases the lA/mcm-IAcm phase 
transition is of the first order. Similar phenomena are ob¬ 
served for the Fmmrn-Frnrn2(ll) phase transition (the 
hysteresis region width is 0.48 kbar for p || [110] and 
0.41 kbar for p || [001]) and for the Fmm2{ll)-Cm phase 
transition (the hysteresis region width is 0.39 kbar for 
p II [110] and 0.62 kbar for p || [001]). 

The most unusual phase transition in Fig. [^ is the 
Cm —>■ Jma2(II) one which is observed in biaxially 


stretched films. As follows from Fig.|^ an abrupt change 
in polarization at this transition is an evidence of the 
first-order phase transition. This conclusion is consistent 
with the absence of any critical phonon mode softening on 
both sides of this transition (Fig. |^. Interestingly, that 
at pressures below the transition pressure, both struc¬ 
tures [Cm and /ma2(II)] satisfy the stability criterion, 
which means that the Jma2(II) phase with a higher en¬ 
thalpy is metastable at these pressures. We attribute the 
transition under discussion to a “sliding off” of the struc¬ 
ture, which occurs when the enthalpy of the Ima2{lV) 
phase becomes lower than that of the Cm phase: thanks 
to the closeness of structures, the height of the potential 
barrier separating them is low. 

It is interesting to also consider the chain of phase tran¬ 
sitions occurring on a branch starting from the /d/mcm 
phase which is metastable at tensile stresses (Fig. [^.1^ 
The calculations of the phonon spectra for the lA/mcm 
phase show that not only the Eu mode is softened with 
increasing tensile stress, but also the Eg mode is soft¬ 
ened, and after the phase transition to the Ima2 phase 
at —4.1 kbar it splits into two components (Fig. [^. At 

— 13.8 kbar, the frequency of one of its components van¬ 
ishes. Adding of the distortions corresponding to this 
unstable mode to the Ima2 structure causes its relax¬ 
ation into a phase which has the same space group, but 
a different octahedral rotations pattern. This is already 
familiar to us the Jma2(II) phase, which was obtained 
earlier upon stretching the film along the lA/mcm —> 
Fmmm —>■ 7^7717712(11) —> Cm —> /777a2(II) chain. As the 
Eg mode originates from the R 25 mode of cubic SrTiOa, 
the cause of the phase transition at —13.8 kbar is the an- 
tiferrodistortive instability, which results in a change of 
the octahedral rotation axis direction from one normal to 
the film plane in the Ima2 phase to one along the polar 
[110] pseudocubic axis in the /77ia2(II) phase. It is un¬ 
usual here that upon the relaxation, the Ima2 structure 
transforms through an intermediate Cc phase into the 
metastable /r?7a2(II) phase, rather than into the ground- 
state structure [the Cm phase according to Fig.j^b)]. 

Similar phenomena are observed in uniaxially stressed 
SrTiOs films at P[ooi] = +19.7 kbar [Fig. j^a)]. The 
Ima2 structure, which becomes unstable above this pres¬ 
sure, also relaxes into the metastable /77ia2(II) phase 
rather than into the ground-state structure [Fm7772(II) 
or Cm]. 

The metastability effects are also characteristic of the 
/r77a2(II) phase at pjooi] = 15-37 kbar and at P[iio] = 

— (11-22) kbar (Fig. [^. 

The behavior of the soft mode frequency in the 
/77ia2(II) phase when changing the in-plane tensile stress 
is also interesting. As follows from Fig.[^ at the isostruc- 
tural /?77a2-/r?7a2(II) phase transition, the majority of 
modes abruptly change their frequencies (dashed lines 
in the figure), whereas the frequency of the soft mode 
does not vanish when approaching to the phase bound¬ 
ary from the Jma2(II) phase. This means that this 
phase transition is of the first order. This conclusion 
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FIG. 5. (Color online) Frequencies of low-energy phonons 
in stable {Mem, M/mcm) and metastable [/ma2, /ma2(II)] 
phases as a function of p || [110] stress. The frequencies in the 
/ma2(II) phase obtained upon releasing the in-plane stress 
are shown by open symbols. 


is also confirmed by the ability to preserve the stability 
of the Ima2(ll) phase when gradually decreasing tensile 
stress (open symbols in Fig. [^. Thus, in the pressure 
range from —4.1 to —13.8 kbar, a hysteresis region oc¬ 
curs on the considered branch of the phase diagram.l^ 
In this region, both Ima2 and /ma2(II) phases satisfy 
the stability criterion, which means, if one takes into 
account the enthalpies of these phases [Fig. [^b)], that 
both phases are metastable. The unusual property of the 
M/mcm —> Ima2 —>■ /TOa2(II) chain of transformations 
is that the phase transition at P[iio] = —4.1 kbar occurs 
from a nonpolar to a metastable polar phase, which trans¬ 
forms into another isostructural metastable polar phase 
with increasing stress (at —13.8 kbar). 

The properties of polar phases coexisting in the hys¬ 
teresis region differ significantly. At p[iio] = —11 kbar, 
the values of spontaneous polarization are 0.128 C/m^ 
in the Ima2 phase, 0.182 C/m^ in the Ima2{]l) phase, 
and 0.096 C/m^ in the ground-state Ftoto 2(II) phase. 
The band gaps calculated in the LDA approximation are 
1.795 eV in the Ftoto 2(II) phase, 1.825 eV in the Ima2 
phase, and 1.861 eV in the Ima2{\l) phase. 


IV. DISCUSSION 

A. Phase diagrams of strained and stressed SrTiOs 

films 

1. Strained SrTiOs films 

To start, we compare the phase diagrams obtained in 
this work for SrTiOa film grown on a cubic substrate 
(the fixed-strain boundary conditions) with the results 


of previous investigations. 

At comp ressive strain, the results of this paper and pre¬ 
vious work^^E^E^ are in good agreement and differ only 
in the strain values at which the lA/mcm-Mcm phase 
transition occurs. This is not surprising if we take into 
account the difference in the sets of material constants 
used in phenomenological models and the difference in 
the calculation techniques and schemes of the pseudopo¬ 
tential construction used in first-principles calculations. 
The difference in the oq values, at which the extrapolated 
squares of the soft-mode frequencies vanish at two sides of 
this phase transition, was found to be only 0.00024 Bohr 
in our calculations, which indicates that the transition is 
close to the second-order one. 

At tensile strain, the results of our and previous cal¬ 
culations are much more different. In the phenomeno¬ 
logical model,!^ in stretched films at all strains the po¬ 
larization is directed along the <100> pseudocubic axis, 
so that our solutions Ima2, Jma2(II), and Cm are ab¬ 
sent. In the phenomenological modeP^it was shown that 
the Ima2 and /ma2(II) solutions appear at low values 
of the 012 coefficient for the {P^Py + Px^z + PyPz) 
term, and the Ima2 phase is stable only in a narrow re¬ 
gion of negative strains. Both these phases disappear 
with increasing 012 (the solution obtained in this case 
agrees with that of Ref. S), and the Cm phase does not 
appear at all. The optimal value of the ai 2 coefficient 
for which the predictions of the phenomenological model 
are closest to the experiment was found in Ref. [22l In 
the phase-field modeP^ at T —>■ 0 the Ima2 phase is 
absent, and it is not clear whether the Cm solution at 
high strains can be associated with our /ma2(II) phase. 
Our results also differ from those of the phenomenolog¬ 
ical model,Pi^ in which its own set of material constants 
was used. According to Ref. HU when the lattice param¬ 
eter ao is increased, the ground-state structures change 
as follows: PAmm -T- 14cm -T I4/mcm -T Ima2 -T 
Fmm2{\l) -T Cm -T /ma2(II) —> Amm2, whereas our 
calculations predict the sequence 14cm —>■ I4/mcm —> 
Ima2 —> Cm -T- Fmm2{\l) -T- /TOa2(II) and that the 
P4mm and Amm2 phases never become the ground-state 
structures. As for the comparison with the results of 
first-principles calculations that took into account both 
order parameters,^^ in the cited work the existence re¬ 
gions of the Ima2 and Jma2(II) phases are separated by 
an intermediate Fmm2{)l) phase. In our calculations, 
the sequence of the ground-state structures for stretched 
films differs from the data of Ref. [18] in that the transi¬ 
tion between the Ima2 and Cm phases occurs abruptly, 
and the Fmm2{ll) phase is the ground-state structure 
only in a narrow region separating the Cm and Ima2{lT) 
phases (Fig. [^. 

2. Stressed SrTiOs films 

We now consider the differences between the phase di¬ 
agrams obtained under the fixed-stress boundary condi- 
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tions for two ways of applying stress and the phase dia¬ 
grams obtained above under the fixed-strain conditions. 

In biaxially compressed film {p || [110]), the transi¬ 
tion pressures of the lAcm-IA/mcm and IA/mcm-Ima2 
phase transitions are close to the values of internal stress 
at the transition points in films grown on a cubic sub¬ 
strate (with a maximum deviation of 0.3-0 .6 kbar). In 
contrast, at the biaxial stretching, the phase diagrams 
differ much stronger [compare Fig. and Fig. Sb)]. 
This is because under tension, the Fmm2(II), Cm, and 
/ma2(II) phases (in which spontaneous strains caused 
by polarization and octahedral rotations lie in the film 
plane) become energetically more favorable. At low 
stresses, a configuration in which the two components of 
spontaneous strain are perpendicular to each other is the 
most favorable, and at high stresses, a configuration in 
which the spontaneous strains are parallel to each other 
and are directed along the [ 110 ] pseudocubic axis is the 
most favorable. Due to an additional contribution to the 
enthalpy proportional to ei -I- 62 = —(where a 
“two-dimensional” elastic compliance is the sum of 
•Sir + S 22 + 25'i 2 elastic compliance tensor components in 
Voigt notation), the phase with the highest is char¬ 

acterized by the strongest decrease of the enthalpy with 
increasing stress, and so it is asymptotically the most 
stable. In our case, this is the /ma2(II) phase. 

Under uniaxial tension {p j| [001]), the phases in which 
the spontaneous strain is directed along the z axis (under 
tension) or those in which the spontaneous strain lies in 
the film plane (under compression) should be more en¬ 
ergetically favorable. However, we should note that the 
condition for the occurrence of a phase transition un¬ 
der uniaxial stress differ significantly from the condition 
for the case of biaxial stress. For example, the lattice 
parameter gq at which the lA/mcm-IAcm phase transi¬ 
tion occurs in SrTiOa is 7.2785 Bohr for p j| [110] and 
7.3018 Bohr for p j| [001]. This reduces the transition 
pressure into the lAcm phase and strongly increases the 
transition pressures into the F'mm2(B) and Cm phases. 
As an additional contribution to the enthalpy for p jj [001] 
is proportional to £3 = —P[ooi]> 5 ' 33 , in this case asymp¬ 
totically the most stable phase is that with the highest 
S 33 value. This phase is also the /ma2(B) phase, but be¬ 
cause the transition pressure into this phase under uniax¬ 
ial compression is sufficiently high (~85 kbar), this phase 
does not become the ground state in the pressure range 
considered in this work [Fig. i:a)]. 

As we have seen, the transition from the fixed-strain 
to the fixed-stress boundary conditions causes a change 
in the order of phase transitions: they become the first- 
order ones, and the width of the corresponding hystere¬ 
sis regions depend on the way of applying stress. We 
attribute the change in the order of phase transitions 
to a well-known renormalization of the coefficients of the 
fourth-order terms in the Landau expansion produced by 
electrostriction, which results in that the second-order 
phase transitions can become the first-order ones. The 
fact that the lattice parameters (as well as the atomic 


positions) near the phase transition points are signifi¬ 
cantly different for two ways of applying stress enables 
to explain why the width of the hysteresis regions near 
these points depends on how the stress is applied. It 
is important to note that the mere fact of appearance 
of the first-order phase transitions in SrTiOa means that 
the phenomenological description of its properties should 
take into account at least the sixth-order terms in the 
power series expansion of the thermodynamic potential. 


B. The nature of metastability effects 

The most interesting result of this work is the observa¬ 
tion of metastable phases and phase transitions between 
them in strained and stressed SrTiOa films. We believe 
that the appearance of these phenomena is due to com¬ 
peting instabilities in this material. 

Indeed, these phenomena were absent in the first- 
principles calculations in which the octahedral rota¬ 
tions were neglected.l^^®^ These calculations predicted 
the second-order PA/mmm —>■ PAmm phase transi¬ 
tion in biaxially compressed films and the second-order 
PA/mmm —>■ Amm2 one in biaxially tensile films. Ac¬ 
cording to our calculations, in which octahedral rota¬ 
tions were neglected, at p = 0 the most stable phase 
is i?3m,^ and the biaxial strain transforms it via an in¬ 
termediate Cm phase, respectively, to the PAmm phase 
(at oq « 7.341 Bohr) or to the Amm2 phase (at oq « 
7.361 Bohr).l^ No signs of metastability were observed 
in these calculations. 

When both ferroelectric and antiferrodistortive insta¬ 
bilities in SrTiOa are taken into account, the metastabil¬ 
ity appears already under the fixed-strain boundary con¬ 
ditions. It is observed in the Fmm2{ll) and Ima2(lT) 
phases at oq = 7.39-7.40 Bohr and in the Pc phase, to 
which an unstable /ma2(B) phase relaxes at Oq = 7.32- 
7.34 Bohr upon condensation of two unstable phonons at 
the center and at the boundary of the Brillouin zone (the 
energy of the Pc phase is significantly higher than that 
of the ground-state lA/mcm or Ima2 phases, Fig. [^. 
The signs of metastability can also be observed in the 
region oq = 7.35-7.37 Bohr where two different solutions 
with the Cm symmetry and different orientation of the 
polarization, which both satisfy the stability criterion, 
appear as a result of relaxation of unstable /ma2(B) and 
Fmm2{ll) phases. 

A particularly large number of metastable phases ap¬ 
pears in SrTiOa under the fixed-stress boundary condi¬ 
tions. Under these conditions, the majority of phases, 
which were simply unstable under fixed-strain condi¬ 
tions, become metastable (in Fig. the points corre¬ 
sponding to these phases are connected by solid lines). 
An analysis shows that the octahedral rotation patterns 
in metastable phases are strongly different from those in 
the ground states. For example, for the Ima2 phase, 
which is metastable under tensile stress, the transition 
to the ground state requires to change the direction of 


the octahedral rotation axis from [001] to [110], and for 
the Pc and Fmmm phases, which are metastable un¬ 
der compressive stress, the transition requires to change 
the direction from [Her] to [001] and from [010] to [001], 
respectively. This suggests that each polar phase has 
its optimal octahedral rotation pattern which stabilizes 
this phase and creates a potential barrier preventing this 
phase to be easily transformed into another structure. 
Such a transformation needs to simultaneously shift a 
large number of the oxygen atoms involved in the octa¬ 
hedral rotations and to change the lattice parameters, so 
that the appearance of metastability is closely related to 
the antiferrodistortive instability. The height of the po¬ 
tential barrier separating the Ima2 and /ma2(ll) phases 
(which can be transformed into each other via an in¬ 
termediate Amm2 phase) is ^24 meV per formula unit 
at P[iio] = “11 kbar. The weakening of the metasta¬ 
bility effects in films grown on cubic substrates, in which 
spontaneous strain cannot be fully realized, confirms this 
explanation. Nevertheless, the competition of the local 
distortions accompanying the appearance of polarization 
and octahedral rotations can be observed even in epitax¬ 
ial films. 

It should be noted that the metastability effects re¬ 
vealed in SrTiOa films are qualitatively different from 
the effects typical for first-order phase transitions. In 
the latter case, the metastability appears in the hystere¬ 
sis region where two phases described by different values 
of the same order parameter fulfill the stability criterion 
and so can coexist. In SrTiOa, the metastability effects 
are due to a complex interaction between the polariza¬ 
tion and the rotational order parameter via a striction 
mechanism. This explains why the obtained solutions 
turn out to be very sensitive to the mechanical boundary 
conditions. 

The consequence of the metastability effects in sys¬ 
tems with competing instabilities is the impossibility, in 
certain situations, to establish whether the obtained so¬ 
lution is the ground-state structure of a system, or just 
a metastable state. As shown in this work, the approach 
used in most first-principles investigations of phase tran¬ 
sitions (according to which the ground state is searched 
by a successive approaching to the structure that sat¬ 
isfies the stability criterion) often results only in the 
metastable states. The ground state in these systems 
can be found only after exploring the entire phase dia¬ 
gram and all transitions between the phases. This anal¬ 
ysis, however, can be facilitated by studying the phase 
diagrams under the fixed-strain boundary conditions. 

It can be seen that the existence of metastable phe¬ 
nomena in SrTiOa results in the appearance of bistability 
regions on the phase diagrams, in which stable and one or 
more metastable polar phases coexist. Let us discuss one 
interesting possibility of using this bistability. The phase 
sequence obtained on the IAcm-I4:/mcm-Ima2 branch 
in stretched SrTiOa films [Fig. [^b)] enables to offer one 
more possible application of such films in non-volatile 
phase change memory devices. Indeed, by applying a 


specific sequence of stresses when cooling the film, a ho¬ 
mogeneous metastable state in the hysteresis region (the 
Ima2 phase) can be prepared in it. Then, the regions 
of equilibrium [FTOm2(II), Cm] or the other metastable 
[Jma2(II)] phases can be formed on the film surface by 
the local optical heating. Due to the bistability of the 
system, both phases can coexist at low temperatures for 
a long enough time, and a contrast in their optical prop¬ 
erties can be used for nondestructive read-out of recorded 
information. Erasing of information can be easily real¬ 
ized by heating the film at zero stress. 

The sequences of the ground states in SrTiOa films 
obtained in this work as a function of applied strain and 
stress correspond to the case of T —)■ 0. At higher tem¬ 
peratures, the phase diagrams may become more com¬ 
plex because of the appearance of domain-like structures. 
Such a structure, which consisted of thick 90° domain 
walls, was revealed in a first-principles-based effective 
Hamiltonian molecular dynamics study of BaTiOa films 
subjected to tensile strain in a wide temperature range 
below the Curie temperature.^^ In unstrained BaTiOa, 
the existence region of this multidomain state coincided 
with that of the orthorhombic phase in bulk BaTiOa. 
However, it remains unclear whether this multidomain 
structure is a real (static) one or just a correlated ther¬ 
mal motion of atoms because the presented data were 
averaged over a very short (4 ps) time interval. The exis¬ 
tence of competing instabilities in SrTiOa can complicate 
its domain structure. Its complex dynamic behavior and 
the metastability effects revealed in this work can result 
in enhanced loss tangent usually observed in ferroelectric 
phases because they both exhibit a retarded reaction to 
an applied electric field. 


V. CONCLUSIONS 

The first-principles calculations of the phase dia¬ 
grams of strained and stressed SrTiOa films have re¬ 
vealed a number of previously unknown metastability 
effects which manifest themselves as coexistence of sev¬ 
eral phases all satisfying the stability criterion. These 
effects become particularly noticeable under the fixed- 
stress boundary conditions. It was shown that the cause 
of these effects is the existence of competing instabilities 
in this material. This suggests that similar phenomena 
may occur in other ferroelectrics with competing insta¬ 
bilities such as NaNbOa and BiFeOa. 

The metastability effects are also important for under¬ 
standing the phenomenon of antiferroelectricity because 
competing instabilities are very typical for antiferroelec- 
tric materials such as PbZrOa. In particular, the electric- 
field-induced polar i?3c phase in PbZrOa is in fact a 
metastable phase, in which the octahedral rotation pat¬ 
tern differs from that in the ground-state nonpolar Pbam 
structure. 

As for the results specific for SrTiOs, it is clear 
that the phenomenological description of its properties 


needs to take into account the higher-order terms in the 
Landau expansion, and so the results of previous studies 
within this approach require a further examination. 
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Appendix: Supplemental material 

The energies and enthalpies of different phases calcu¬ 
lated for a given strain or stress are compared in Ta¬ 
bles A1 and A2. The optimized lattice parameters and 
atomic positions for all structures considered in this work 
are given in Tables A3 to A12. 
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TABLE Al. Energies of different phases in strained SrTiOa film fixed on a square substrate with a lattice parameter ao (in 
meV). The energy of high-symmetry PA/mmm phase is taken as the energy reference. The energies of stable phases are in 
boldface. 


ao 

(Bohr) 

lA/mcm 

74cm 

Fmmm 

Ima2 

Phases 

Fmm2 

7ma2(II) 

Fmm2{ll) 

Cm 

7.16 

-93.03 

- 95.72 

— 

— 

— 

— 

— 

— 

7.20 

-73.46 

- 74.53 

— 

— 

— 

— 

— 

— 

7.22 

-64.97 

- 65.53 

— 

— 

— 

— 

— 

— 

7.24 

-57.19 

- 57.47 

— 

— 

— 

— 

— 

— 

7.26 

-50.17 

- 50.28 

— 

— 

— 

— 

— 

— 

7.27 

-46.95 

- 46.97 

— 

— 

— 

— 

— 

— 

7.275 

-45.40 

- 45.43 

— 

— 

— 

— 

— 

— 

7.28 

- 43.91 

— 

— 

— 

— 

— 

— 

— 

7.29 

- 41.01 

— 

— 

— 

— 

— 

— 

— 

7.30 

- 38.28 

— 

— 

— 

— 

— 

— 

— 

7.31 

- 35.70 

— 

— 

— 

— 

— 

— 

— 

7.315 

- 34.45 

— 

— 

— 

— 

— 

— 

— 

7.32 

-33.24 

— 

-26.67 

- 33.26 

-33.23 

-27.67 

— 

-28.21 

7.33 

-30.94 

— 

-26.89 

- 31.07 

-30.99 

-27.81 

— 

-28.03 

7.34 

-28.76 

— 

-27.12 

- 29.09 

-28.96 

-27.96 

-27.40 

-28.00 

7.345 

-27.71 

— 

-27.24 

- 28.19 

-28.01 

-28.03 

-27.63 

-28.05 

7.35 

-26.69 

— 

-27.36 

-27.37 

-27.12 

- 28.17 

-27.87 

- 28.17 

7.36 

-24.71 

— 

-27.61 

-25.89 

-25.47 

-28.48 

-28.39 

- 28.49 

7.37 

-22.86 

— 

-27.86 

-24.65 

-24.04 

-28.97 

-28.99 

- 29.04 

7.38 

-20.98 

— 

-28.02 

-23.49 

-22.65 

-29.44 

-29.51 

- 29.54 

7.39 

-19.38 

— 

-28.30 

-22.77 

-21.62 

-30.25 

- 30.33 

— 

7.40 

-17.87 

— 

-28.58 

-22.22 

-20.77 

-31.21 

- 31.21 

— 

7.42 

-15.05 

— 

-29.07 

-21.91 

-19.61 

- 33.55 

-33.19 

— 

7.44 

-12.59 

— 

-29.56 

-22.53 

-19.29 

- 36.65 

-35.56 

— 
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TABLE A2. Entahlpies of different phases in [110]-stressed SrTiOa film (in meV). The energy of high-symmetry PA/mmm 
phase is taken as the energy reference. The enthalpies of stable phases are in boldface. 


P[110] 

(kbar) 

/4cm 

/4/mcm 

Ima2 

/ma2(II) 

Phases 

Imma 

Fmmm 

Fmm2{ll) 

Cm 

44.13 

— 

— 

-23.08 

-40.48 

— 

-32.40 

-36.86 

— 

36.78 

— 

— 

-22.09 

-36.67 

— 

-32.15 

-34.88 

— 

29.42 

— 

-19.36 

-22.21 

-33.86 

— 

-31.96 

-33.42 

-33.80 

22.07 

— 

-21.95 

-23.26 

-31.79 

— 

-31.70 

-32.33 

-32.41 

18.39 

— 

-23.28 

-24.11 

-30.98 

— 

-31.53 

-31.90 

-31.91 

14.71 

— 

-24.72 

-25.15 

-30.34 

— 

-31.39 

-31.5eF| 

-31.54 

11.03 

— 

-26.30 

-26.48 

-29.95 

— 

-31.39 

-31.45 

— 

7.36 

— 

-27.76 

-27.81 

-29.53 

-29.11 

-31.21 

-31.24 

— 

5.52 

— 

-28.52 

-28.54 

-29.37 

-29.07 

-31.10 

-31.12 

— 

4.60 

— 

-28.91 

-28.91 

-29.28 

-29.05 

-31.06 

-31.08 

— 

3.68 

— 

-29.30 

— 

-29.22 

-29.02 

-31.02 

— 

— 

0.0 

— 

-30.88 

— 

-28.97 

-28.92 

-30.88 

— 

— 

-3.68 

— 

-32.45 

— 

-28.80 

— 

-30.81 

— 

— 

-7.36 

— 

-34.12 

— 

-28.70 

— 

-30.67 

— 

— 

-11.03 

— 

-35.86 

— 

-28.59 

— 

-30.47 

— 

— 

-14.71 

-37.69 

-37.63 

— 

-28.50 

— 

-30.35 

— 

— 


^ Response function calculations indicate an instability of this phase against its transformation into the Cm phase. 
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TABLE A3. Lattice parameters, Wyckoff positions, and 
atomic coordinates for the I^jmcm phase of SrTiOa at p = 0 
(the rotations are around the [001] pseudocubic axis). 


a = 10.333105 Bohr, c = 14.803747 Bohr 


Atom 

WP 

X 

y 

Z 

Sr 

4b 

0.000000 

0.500000 

0.250000 

Ti 

4c 

0.000000 

0.000000 

0.000000 

O 

4a 

0.000000 

0.000000 

0.250000 

O 

Sh 

0.214346 

0.714346 

0.000000 


TABLE A4. Lattice parameters, Wyckoff positions, and 
atomic coordinates for the 14cm phase of SrTiOa film fixed 
on a substrate with ao = 7.20 Bohr (the rotations are around 
the [001] pseudocubic axis, the polarization is along the same 
axis). 


a = 10.182338 Bohr, c = 14.996956 Bohr 


Atom 

WP 

X 

V 

z 

Sr 

4b 

0.500000 

0.000000 

0.253951 

Ti 

4a 

0.000000 

0.000000 

0.008361 

O 

4a 

0.000000 

0.000000 

0.249368 

O 

8c 

0.207733 

0.707733 

-0.000840 

TABLE 

A5. Lattice parameters. 

Wyckoff positions, and 


atomic coordinates for the Fmmm phase of SrTiOa film fixed 
on a substrate with ao = 7.38 Bohr (the rotations are around 
the [010] axis). 


a = 14.760000 Bohr, c = 14.593072 Bohr 


Atom 

WP 

X 

y 

z 

Sr 

8h 

0.000000 

0.250507 

0.000000 

Ti 

8d 

0.250000 

0.000000 

0.250000 

O 

8i 

0.000000 

0.000000 

0.216530 

O 

8/ 

0.250000 

0.250000 

0.250000 

O 


0.283971 

0.000000 

0.000000 


TABLE A6. Lattice parameters, Wyckoff positions, and 
atomic coordinates for the Fmm2 phase of SrTiOa film fixed 
on a substrate with ao = 7.38 Bohr (the polarization is along 
the [100] pseudocubic axis, the rotations are around the [001] 
axis; the F2mm setting). 



a = 14.760000 Bohr, c = 

14.703578 Bohr 


Atom 

WP 

X 

V 

Z 

Sr 

8d 

0.006484 

0.000000 

0.249962 

Ti 

8c 

0.258024 

0.249145 

0.000000 

O 

8c 

-0.000988 

0.281480 

0.000000 

O 

4a 

0.218214 

0.000000 

0.000000 

O 

4a 

-0.219248 

0.000000 

0.000000 

O 

8b 

0.246996 

0.250000 

0.250000 


TABLE A7. Lattice parameters, Wyckoff positions, and 
atomic coordinates for the Fmm2{ll) phase of SrTiOs film 
fixed on a substrate with ao ~ 7.38 Bohr (the polarization is 
along the [100] pseudocubic axis, the rotations are around the 
[010] axis; the F2mm setting). 



a = 14.760000 Bohr, c = 

14.591528 Bohr 


Atom 

WP 

X 

y 

z 

Sr 

8c 

0.006027 

0.250548 

0.000000 

Ti 

8d 

0.258045 

0.000000 

0.250988 

O 

8d 

0.499078 

0.000000 

0.283937 

O 

8b 

0.247904 

0.250000 

0.250000 

o 

4a 

0.283166 

0.000000 

0.000000 

o 

4a 

-0.285274 

0.000000 

0.000000 


TABLE A8. Lattice parameters, Wyckoff positions, and 
atomic coordinates for the Ima2 phase of SrTiOs film fixed 
on a substrate with ao = 7.38 Bohr (the polarization is along 
the [110] pseudocubic axis, the rotations are around the [001] 
axis; the I2am setting). 



a = 10.436896 Bohr, c = 

14.702085 Bohr 


Atom 

WP 

X 

y 

z 

Sr 

4a 

0.511582 

0.000000 

0.000000 

Ti 

4b 

0.014611 

0.001526 

0.250000 

O 

4a 

-0.004529 

0.000000 

0.000000 

O 

4b 

0.717532 

0.218809 

0.250000 

o 

4b 

-0.219196 

0.717984 

0.250000 


TABLE A9. Lattice parameters, Wyckoff positions, and 
atomic coordinates for the /ma2(II) phase of SrTiOs film 
fixed on a substrate with ao = 7.38 Bohr (the polarization 
is along the [110] pseudocubic axis, the rotations are around 
the same axis; the 12am setting). 



a = 10.436896 Bohr, c = 

14.594969 Bohr 


Atom 

WP 

X 

y 

Z 

Sr 

4b 

0.510431 

0.001862 

0.250000 

Ti 

4a 

0.012530 

0.000000 

0.000000 

O 

4b 

0.497252 

0.547074 

0.250000 

O 

8c 

0.250394 

0.249830 

0.024560 
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TABLE AlO. Lattice parameters, Wyckoff positions, and 
atomic coordinates for the Cm phase of SrTiOa film fixed 
on a substrate with ag = 7.36 Bohr (the polarization and the 
rotation axis are in the xy plane; the Am setting). 



a — 10.408612 Bohr, c = 

14.619606 BohT 

1 

Atom 

WP 

X 

y 

z 

Sr 

2a 

0.504500 

0.503322 

0.000000 

Sr 

2a 

0.505691 

-0.000127 

0.000000 

Ti 

46 

0.005456 

0.002834 

0.250500 

O 

2a 

0.029421 

0.454223 

0.000000 

O 

2a 

-0.029241 

0.043574 

0.000000 

O 

46 

0.001801 

0.248585 

0.280345 

O 

46 

0.501557 

0.248085 

0.734216 


The translation vectors are [100], [110], and [001]. 


TABLE All. Lattice parameters, Wyckoff positions, and 
atomic coordinates for the Fmm2(III) phase of SrTiOs film 
subjected to a stress of P[ooii = —14.7 kbar (the polarization 
is along the [001] pseudocubic axis, the rotations are around 
the [010] axis). 


a = 14.59641 Bohr, 6 = 14.77120 Bohr, 

c = 14.72473 Bohr 

Atom 

WP 

X 

y 

Z 

Sr 

8c 

0.000000 

0.249575 

0.001152 

Ti 

8d 

0.249152 

0.000000 

0.252941 

O 

4a 

0.000000 

0.000000 

0.210624 

O 

4a 

0.000000 

0.000000 

0.279767 

O 

86 

0.250000 

0.250000 

0.244342 

O 

8d 

0.284337 

0.000000 

-0.004654 


TABLE A12. Lattice parameters, Wyckoff positions, and 
atomic coordinates for the Pc phase of [001]-stressed SrTiOs 
film (p[ooi] = —14.7 kbar); the Pa setting. 


a = 10.351368 Bohr, b = 10.405403 Bohr, 



C — 

14.747770 Bohr, 

j3 = 89.9952° 


Atom 

WP 

X 

y 

Z 

Sr 

2a 

0.507791 

0.252513 

0.255092 

Sr 

2a 

0.497563 

0.247451 

0.755063 

Ti 

2a 

0.002976 

0.248822 

0.007148 

Ti 

2a 

0.002331 

0.251203 

0.507150 

O 

2a 

0.001328 

0.204662 

0.249255 

O 

2a 

0.004428 

0.295355 

0.749257 

O 

2a 

0.244962 

-0.007996 

-0.024935 

O 

2a 

0.760337 

0.507554 

0.475052 

O 

2a 

-0.738300 

0.491565 

0.022296 

O 

2a 

-0.255419 

0.007896 

0.522307 
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